TRANSCENDENTAL HARMONIC MAPPINGS AND 
GRAVITATIONAL LENSING BY ISOTHERMAL 

GALAXIES 



DMITRY KHAVINSON AND ERIK LUNDBERG 



Abstract. Using the Schwarz function of an ellipse, it was re- 
cently shown that galaxies with density constant on confocal el- 
hpses can produce at most four "bright" images of a single source. 
The more physically interesting example of an isothermal galaxy 
has density that is constant on homothetic ellipses. In that case 
bright images can be seen to correspond to zeros of a certain tran- 
scendental harmonic mapping. We use complex dynamics to give 
an upper bound on the total number of such zeros. 



1. Introduction 

In this paper, we obtain an upper bound for the number of solutions 
of the equation 



'1.1) arcsin 



z + w ^ 

where w is a complex parameter, and A; is a real parameter. 



Our motivation for doing so is that solutions of (1.1) in fact corre- 
spond to virtual images observed when the light from a distant source 
passes near an isothermal, ellipsoidal galaxy. Indeed, using the complex 
formulation of the thin-lens approximation ([IH]), the lensing equation 
is calculated by finding the Cauchy transform of the mass distribution 
projected to the "lens plane" . This was carried out in [7J with the fol- 
lowing result (we also sketch the derivation in the last section for the 
reader's convenience). 



^1.2) C arcsin [ ^ ] + ui = ( 



^Both authors gratefully acknowledge partial support from the National Science 
Foundation. 
Date: 2009. 
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Here, we take the principal branch of arcsin, C and c are real con- 
stants depending on the elliptical projection of the galaxy onto the lens 
plane, and u is the position of the source (projected to the lens plane). 



and values of ( which satisfy (1.2) give positions of the observed im- 
ages. Changing variables to z = w = uj/C, and k = c/C puts 



(1.2) into the form of equation (1.1) while preserving the number of 



solutions. 

We should mention that the anti-analytic potential in the lensing 
equation considered here (and, also, in p| and [7]) differs from the po- 
tential in the lensing equation in the model often used by astrophysi- 
cists (see |8] and the references therein), where the projected mass 
density is supported in the entire complex plane. Both models use the 
"isothermal" density proportional to 1/t on ellipses {x^ / a^+y"^ jlP' = t^} 



(a and b fixed). The model considered here that yields equation (1.2) 
assumes that the density is zero for all t greater than some value (see 
Appendix). Letting the density have infinite support assumes that the 
galaxy has infinite mass and fills the universe, yet it is the simplest way 
to avoid giving the galaxy a "sharp edge" and astronomers have found 
that the model behaves reasonably in the region where the lensed im- 
ages occur. We consider the model with physically realistic compact 
support but less realistic "sharp edge" for a mathematical reason: in 



that setting, lensed images described by solutions of equation (1.1) 
correspond to zeros of a harmonic function (We note that models with 
"sharp edges" have been considered by astrophysicists as well, cf. the 
recent preprints [13] and P3]). 

For gravitational lenses consisting of n point masses, Mao, Fetters, 
and Witt [I3] suggested (1997) that the bound for the number of im- 
ages was linear in n (Bezout's theorem provides a bound quadratic in 
n). Rhie refined this in 2001, conjecturing that a gravitational lens 
consisting of n point masses cannot create more than 5n — 5 images 
of a given source [16]. In 2003, she constructed point-mass configura- 
tions for which these bounds are attained [T7]. The first author and 
G. Neumann [TU] settled her conjecture by giving a bound of 5n — 5 
zeros for harmonic mappings of the form r{z) — z, where r{z) is ra- 
tional of deg n > 1. (See [TT] for the exposition and further details.) 



Solutions of (1.1) are zeros of a transcendental harmonic function, so 
extending the techniques used in [10] will require some care (a priori, 
it is not even clear that the number of zeros is finite, cf. [1], [12]). Still, 
our approach draws on the same two main results: (i) the argument 
principle generalized to harmonic functions and (ii) the Fatou theorem 
from complex dynamics regarding the attraction of critical points. In 
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the next section, we will formulate (i). (ii) will have to be modified for 
our purposes, so ideas from complex dynamics are worked from scratch 
into the proof of Lemma |3.3| in the third section. 

Acknowledgement: We would like to thank Walter Bergweiler, 
Alex Eremenko, and Charles R. Keeton for stimulating discussions and 
for sharing some of their unpublished work with us. 

2. Preliminaries: The Argument Principle 

In order to state the generalized argument principle (see [5] for a 
complete exposition and proof), we need to define the order of a zero 
or pole of a harmonic function. A harmonic function h = f + g, where 
/ and g are analytic functions, is called sense-preserving at zq if the 
Jacobian Jh{z) = — > for every z in some punctured 

neighborhood of Zq. We also say that h is sense-reversing if h is sense- 
preserving at Zq. If h is neither sense-preserving nor sense- reversing 
at Zq, then zq is called singular and necessarily (but not sufficiently) 
Jh{zQ) = 0, cf. [5], Ch. 2. The order of a non-singular zero is given by 
2^Ac arg /i(z), where C is a sufficiently small circle around the zero. 
The order is positive if h is sense-preserving at the zero and negative 
if h is sense-reversing. 

Suppose h is harmonic in a punctured neighborhood of Zq. We will 
refer to pole of h if h{z) oo as z ^ zq. Let C be an oriented 

closed curve that contains neither zeros nor poles of h. The notation 
2^ Ac arg h{z) denotes the increment in the argument of h{z) along C. 
Following [T9], the order of a pole of h is given by —i^AcS3:gh{z), 
where C is a sufficiently small circle around the pole. We note that if 
h is sense-reversing on a sufficiently small circle around the pole, then 
the order of the pole will be negative. We will use the following version 
of the argument principle which is taken from [I^ : 

Theorem 2.1. Let F be harmonic, except for a finite number of poles, 
in a Jordon domain D. Let C be a curve contained in D not passing 
through a pole or a zero, and let R be the open, bounded region sur- 
rounded by C. Suppose F has no singular zeros in R and let N be the 
sum of the orders of the zeros of F in R. Let P be the sum of the 
orders of the poles of F in R. Then Ac a.igF{z) = 2tt{N — P). 

3. An Upper Bound for the Number of Images 



Lemma 3.1. The solutions of equation (1.1) are all contained in a 
rectangle, R := {\Re{z)\ < 7i/2, \Im{z)\ < M}, where M is sufficiently 
large. 
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Proof. The requirement that |i?e(2;)| < 7r/2 is immediate since this strip 
is the image of C under the principal branch of arcsin. To see that there 
exists an M such that solutions of satisfy \Im(z)\ < M, take sin 
of both sides. This leads to 



(3.1) =^=^ + ^. 

sm[z) 

We consider the modulus of each side of (3.1) ioi z = x + iy with 
large values of \y\. Recall, sm{x+iy) = sin(x) cosh{y)+i cos(x) sinh(?/). 
As y ±oo, I sm{x+iy) I ~ A;/ sin^ x cosh^ y + cos^ x sinh^ y 0, imi- 
formly in x. On the other hand, |;s + w| oo. 

□ 

Remark: By this lemma, we can bound the number of solutions of 



(1.1 ) by bounding the number of zeros of F(z) := z + w — , , in the 



rectangle, R. Let us calculate the increment of the argument of F(^z) 
when dR is traced counterclockwise. 

Lemma 3.2. Aa^argF(2;) > ~2ix , where F{z) := z \ w — 



sin(z) 



Proof. Consider the four links 14(t) = ±f ± i{-M + t), < t < 2M, 
Pl±{t) = ±(| — t) ±iM, < t < TT, which trace the right, left, top, and 
bottom edges, respectively. We need to determine the effect of the term 
— ^=^. Without this term, F(z) is iust translation z ^ z + w, and 

Sin(2) / \ y J 

in that case F{V±{t)) and F{H±{t)) trace the edges of the translated 
rectangle. 

By choosing M large enough in the previous lemma, we can neglect 
the term — ==^ on the top and bottom edges. On the right edge, 

sin(z) ^ ^ I 

„^ , = — , ,x is pure real and increases monotonically from 

sin(V+(t)) cosh(-M+t) J" J 

a small value at t = to the value k aX t = M. On the interval 
M < t < 2M, eosh(-^M+t) decreases monotonically from k aX t = M 
back to the original value at t = 2M. Similarly, on the left edge, 

> —k. Thus, the effect of the term — . ^, , is to 



sin(y_(i)) cosh(-M+t) ' ' 

bend the left and right sides of the translated rectangle inward, so that 
they cross each other if and only if A; > | (compare the two images in 
figure [T]) . 

If k < 71/2, then the images of the left and right edges do not inter- 
sect, and either Ag^aig F{z) = 27i, or, if F{dR) does not surround the 
origin, AgjiaTgF^z) = 0. See the left image in figure [l] 
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Figure 1. The image of OR under F{z) with w = and 
two choices for the value oi k. In the left, = 1 < 7r/2. 
In this case, Aa_RargF(z) = 1. If we set w to, say, 1 then 
we have AQRaigF{z) = 0. In the right, k = 2 > 7r/2, 
and AdRa.YgF{z) = —1. If we set w to, say, 1 or i then 
we have AQ^aigF^z) = or 1, respectively. 

If k > tt/2, then the images of the left and right edges intersect ex- 
actly twice. In this case, there is a third possibility in which Ag/j arg F{z) 
= —2n. See the right image in figure [TJ 

□ 

Define the anti-analytic function f{z) := '^=^~w whose fixed points 

coincide with the zeros of F{z). Notice that sense-preserving zeros of 
F{z) coincide with attracting fixed points of f{z) (see [B] and [1]). In- 
deed, a fixed point zq is attracting when |/'(2;o)| < 1 which, if satisfied, 
holds in a neighborhood of zq so that 1 — > near zo which is 

the condition required for F{z) to be sense-preserving at zq. We use 
complex dynamics to bound the number of attracting fixed points of 
f{z). The version of the Fatou theorem found in most textbooks on 
complex dynamics such as j4j falls short, since / has infinitely many 
critical points, and / is not analytic but rather anti-analytic. The 
extensions of the Fatou theorem discussed in, e.g., the survey [2] are 
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almost sufficient to cover our situation, but / is still assumed to be 
analytic. For the reader's convenience, we give an independent direct 
proof of the next lemma. 

Lemma 3.3. The number of sense-preserving zeros, n+, is at most 3. 

Proof. Suppose zq is a sense preserving zero. Then |/'(2;o)| = 7 < 1, 
and there is a neighborhood, Ui, of which is contracted to zq under 
iterations of f{z). Consider /^^(C) = 9'rcsin(==^). Recall that arcsin 
has a single-valued, injective branch in any simply connected domain 
in C \ {±1, oo}. Since =j injectively maps C \ {—w, —w ± k} into 

C \ {±1, oo}, we conclude f~^ has a single-valued, injective branch in 
any simply connected domain which omits {—w, —w ± k}. 

Now we consider the algorithm described in the proof of Fatou's the- 
orem. Suppose Ui omits {—w, —w±k}. Then choose the single-valued 
branch of f~^ so that f~^{zo) = zq. Then U2 '■= /"^(f/i) is a simply 
connected domain, since f~^ is analytic and injective. Therefore, if U2 
omits {—w, —w ± k} then f~^ has a single-valued extension to 1/2- We 
can proceed inductively, defining t/„+i := f^^{Un), provided Un omits 
{—w, —w ± k}. 

If the algorithm does not stop, then it produces an infinite family, 
{/~^"}^Q, of functions analytic on the domain Ui (odd iterates of f~^ 
are anti-analytic). This family omits {—w, —w ± k} and is therefore a 
normal family (by Montel's theorem). This contradicts the divergence 
of ^/~"(C)lc=2o which can be seen by repeated application of the chain 
rule along with the fact that ;^/~^(C)|f=2o has modulus I/7 > 1. Thus, 
at least one of the three points —w, —w ± A; is attracted to zq under 
iteration of f{z), and the lemma follows. □ 

Theorem 3.4. The number of solutions to (1.1) is bounded by 8. 



Proof. By the remark following Lemma 3.1, the total number of so 



lution to (1.1) equals the total number of zeros of F[z) in R. Recall 
that F{z) is called "regular" if it is free of singular zeros (see [9] and 
[6]). Suppose for the moment that F{z) is regular. Then, the total 
number of zeros of F{z) in i? is = n+ + , where n_|_ and n_ count, 
respectively, the sense-preserving and sense-reversing zeros of F{z) in 



R. By Lemma 3.2 and Theorem 2.1 —1 < N — P. By Lemma 3.3 



and the fact that F{z) has one sense-reversing pole in i?, this becomes 
— 1 < 3 — + 1, so that n_ < 5. Thus, = + < 8. 

Fix k. There is a dense set of parameters w for which F{z) is regular. 

I — { — - — ) I = 1 1 under z — ^ 

I dz V sin(2;) ' I J sin(2) ' 



Indeed, consider the image of {z : \4-{^^i-^)\ = 1| under z — ==. This 
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set has empty interior, and if w is in its complement, F{z) is free of 
singular zeros. 

Now suppose F{z) is not regular. If F{z) is perturbed by a suffi- 
ciently small constant, the number of zeros is not changed. By the 
preceding, we can choose an arbitrarily small constant so that the per- 
turbation is regular. Thus, the bound, N < 8, holds for all F{z). □ 



4. Concluding Remarks 




Figure 2. Four images of a light source behind an ellip- 
tical galaxy. (Credit: NASA, Kavan Ratnatunga, Johns 
Hopkins Univ.) 



So far, astronomers have observed only up to 5 images (4 bright + 
1 dim) produced by an elliptical lens (see figure [2]). In |8] there have 
been constructed explicit models (depending on the semiaxes of the 
ellipse) having 9 images (8 bright -|- 1 dim) but only in the presence of 
a shear, i.e. a (linear) gravitational pull from infinity (a term added 



to equation (1.1)). So far, we have not been able to obtain a universal 



bound in the presence of a shear that is similar to Theorem 3^ It 
seemed, based on NASA observations, natural to conjecture that, in 
the absence of shear, there can be at most 4 bright images. Yet, recently 
W. Bergweiler and A. Eremenko generated an example with 6 bright 
images [3]. With their kind permission, we include their example (see 
figure |3). Based on their breakthrough and the investigation carried 
out in [8J, we conjecture the following. 

Conjecture 4.1. (i) The number of bright images lensed by an isother- 
mal elliptical galaxy without shear is at most 6. (ii) In the presence of 
a shear, the number of bright images is at most 8. 

We caution the reader that in [8] the mass density was assumed to 
be extended all the way to infinity, so the lensing potential in [8] was 
different from the one we consider here (and in [3]). 
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Figure 3. Equation has 6 solutions when k=1.92 
and w=-.67i. Choosing a = 1, b = .041 , and M = 2 
(see Appendix) leads to k = 1.92 and gives the picture 
of the six images shown here in the (-plane along with 
the galaxy's elliptical silhouette and the source (plotted 
as box).) 

Added in proof: After this paper was submitted, Walter Berg- 
weiller and Alex Eremenko succeded in proving part (i) of the conjec- 
ture. 

5. Appendix: derivation of the complex lensing equation 
for the isothermal elliptical galaxy 

Suppose that light from a distant source star is distorted as it passes 
by an intermediate, continuous distribution of mass which does not 
deviate too far from being contained in a common plane (the "lens 
plane" ) perpendicular to our hne of sight. Let /i(z) denote the projected 
mass density. Then basic results from General Relativity combined 
with Geometric Optics (see [IB]) lead to the following lensing equation 
relating the position of the source (projected to the lens plane) w to 
the positions of lensed images z. 



(5.1) 
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Consider, first, the case when the projected density fi{z) = D is 

2 2 

constant and supported on := + ^ < 1, a > b > 0}, an elhpse. 
Then equation (5.1) becomes 



DdAjQ 



+ w. 



By the complex Green's formula, for z outside Q (i.e., for "bright" 
images), this becomes 



D 



+ w. 



2^ J an C 

The Schwarz function (by definition, analytic and = ^ on dQ) for 
the ellipse equals (c^ = — 6^): 



2ab 



< 



c 

2ab 



. C + ^(C-v/C^) 

= 5i(C) + ^2(C) 

where Si is analytic in Q, and 5*2 is analytic outside Q and 52(00) = 0. 
Since z is outside fl, combining this with Cauchy's formula gives 

2ab 



SiOdC 



+ w 



rr—D{-z 



c^j + w 



2? Jon C-z 
for the right-hand-side of the lensing equation. 

Next consider the case of "isothermal" density supported on f2, 
M/t on dVtt, ■■= tVt = {"^ + ^ < t^},t < 1, and M a constant. 



Then the Cauchy potential term in the lensing equation (5.1) be- 
comes 



(5.2) 



dAiO 



M 


' d 




L ~t 


jit 





dt 



dAiO 



For the inside integral, we see that L . - — v - — v ^ -,. 

° ' Jilt t(—z JU C—z/t 

which according to our previous calculation is Co{z—^/z^ — cH^), where 
the constant Cq depe nds only on Q. Now the t-derivative of this is 



Co-7p=f^. Thus (5.2) becomes MCq 



dt 
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Finally, we arrive at (1.2), the lensing equation for the isothermal 
elliptical galaxy. 



where C = ^M. 



z = C arcsin (^-^ + w, 
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